Correction to the Chiral Magnetic Effect from axial-vector interaction 
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The recent lattice calculation at finite axial chemical potential suggests that the induced current 
density of the chiral magnetic effect (CME) is somehow suppressed comparing with the standard 
analytical formula. We show in a NJL-type model of QCD that such a suppression is a natural 
result when considering the influence of the attractive axial- vector interaction. We point out that 
the lattice result doesn't need to be quantitatively consistent with the analytical formula due to 
the chirality density-density correlation. We also investigate the nonperturbative effect of instanton 
molecules on the CME. Since an unconventional repulsive axial-vector interaction is induced, the 
CME will be enhanced significantly by the instanton-anti-instanton pairings. Such a prediction needs 
to be tested by more improved lattice simulations. We further demonstrate that the axial-vector 
interaction plays an important role on the T — fiA phase diagram. 



I. INTRODUCTION 



The topological charge of gauge field theory is defined 



as 



Qt = 



32tt 2 



d 4 x Tt[FF], 



(1) 



where F and F refer to the field strength tensor and 
its dual, respectively. It is well-known that Quantum 
Chromodynamics(QCD) contains nontrivial gauge con- 
figurations carrying the topological charge with an inte- 
ger number p|. These configurations interpolate between 
the topologically different vacua of QCD characterized 
by different Chern-Simons numbers. For a long time the 
experimental evidence for the existence of the topologi- 
cal gluon configurations only comes indirectly from the 
meson spectrum[2]. 

Recently, the study of the chiral magnetic 
effect(CME)H Q has received much attentions since 
the topological configurations of QCD may be observed 
directly in the noncentral heavy-ion collisions. The 
essence of the CME is the chiral imbalance induced 
by the nonzero topological charge through the axial 
anomaly of QCD 



(N R - N L ) t -_ 



(N R - N L ) 



L)t=-c 



-2Qi 



(2) 



where Nr(Nl) denotes the quark number of the right- 
handed (left-handed) chirality in the chiral limit. Eq. @ 
indicates the chiral asymmetry is generated by the 
nonzero topological excitations. In a strong magnetic 
field which may be produced in the noncentral heavy- 
ion collisions, a net electric current will be induced along 
the direction of the magnetic field due to the chirality 
imbalance. 

In order to study the effect of the chiral asymmetry, 
usually the axial chemical potential ha is introduced 
in a grand canonical ensemble [4|. The physical mean- 
ing of ha is the difference of the right- and left-handed 



quark chemical potentials. This quantity can be identi- 
fied as the time derivative of the angle of QCD, namely 
[i A = do6/2Nf\^. In the presence of an external mag- 
netic field B, the analytical expression of the induced 
electric current density for the CME at finite ha has been 
given in Q which takes the form 



3t 



CN c N f e 2 ii A B, 



(3) 



where C = \/2n 2 . Note that all quarks carrying the 
same electric charge is assumed in <j3j) . Obviously, this 
equation explicitly breaks the P- and CP-symmetry. 

Equation. (J3]) indicates that the hot quark-gluon 
plasma(hQGP) is an ideal place to study the physics of 
P- and CP-odd excitations of QCD. First, the produc- 
tion probability of the configurations with nonzero wind- 
ing numbers becomes higher since jumping over the po- 
tential barrier between distinct classical vacua is proba- 
ble at high temperatures. The configurations responsible 
for such thermal transitions are called sphaleronsJH, H[ 
in QCD, which might occur at a copious rate at high 
temperature compared to the low-temperature tunnel- 
ing by the instanton. Second, a strong magnetic field 
can be produced in non-central heavy- ion collisions d,0|. 
The study based on the UrQMD modelQ suggests that 
the eB produced at RHIC and LHC can be as large as 
eB w 2m 2 and eB » 15m 2 ,, respectively 1 . So in a non- 
central heavy-ion collision, the induced electric current 
may lead to an excess of the positive electric charge on 
one side of the reaction plane and the negative electric 
charge on the other. Such a charge separation effect may 
be observed experimentally. Recently, a conclusive obser- 
vation of charge azimuthal correlation has been presented 
by the Star collaboration [8j which may result from the 
CME with local V- and CP-violation. 

Theoretically, the investigation of the CME has been 
performed recently in the Lattice QCD simulation by in- 
troducing a finite ha |2| . The advantage of the lattice sim- 
ulation at finite axial chemical potential is that there is 
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1 The magnetic field strength can be translated into the CGS sys- 
tem with the identity m% ~ 10 ls Gauss, where m w = 140 MeV. 
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no sign problem comparing with the case at hnite baryon 
chemical potential. In this study, the CME described 
by Eq.([3]) has been confirmed qualitatively. But at the 
quantitative level, the obtained electric current density 
is significantly suppressed [9[ compared to the analytical 
formula ((3]). The further study at the quenched level of 
QCD suggests that the lattice data are sensitive to the 
lattice spacing[l(| and the obtained data increase signif- 
icantly compared to the ones inQ. However, the devia- 
tion of the new lattice data from the analytical formula is 
still sizable even some systematic errors have been taking 
into account [Toj . 

The suppression of the lattice data may stem from 
some errors of the lattice simulation and/or the non- 
perturbative effects of QCD. Accordingly, besides tak- 
ing more improved lattice calculations, physically under- 
standing the current lattice result through other methods 
of QCD is very important and necessary. On the other 
hand, due to the limitation of the lattice calculation, in- 
vestigating the non-perturbative corrections and making 
new predictions on the CME are also very required by 
using the effective theories or models of QCD. 

The purpose of this paper has two aspects: First, we 
will try to disclose the possible mechanism for the devia- 
tion of the present lattice data from the analytical result 
by taking into account some QCD corrections; Second, 
we will explore the non-perturbative influence of the in- 
stanton molecules on the CME and make some predic- 
tions. Note that both aspects are closely related to the 
effective four-quark interaction in the axial-vector chan- 
nel, where the chirality density-density correlation plays 
an important role. 

The paper is organized as follows: In Section[Hl we give 
the effective quark interactions of QCD for T > T c . In 
section IIII| we present the correction to the CME from 
the axial-vector interaction. Section IIVI is devoted to 
the numerical results and related discussions. The last 
section is the summary and conclusion. 



II. EFFECTIVE QUARK INTERACTIONS FOR 

T>T C 

The experimental data of RHIC suggests that strongly 
coupled quark-gluon plasma(sQGP) is formed for T > T c , 
which exhibits the behavior of a perfect fluid. Now it 
is widely believed that the non-perturbative QCD still 
plays very important roles on the hQGP up to 2 — 3T C . 
Since the chiral magnetic effect may be observed directly 
in heavy-ion collisions, one naturally expects that the 
role of the non-perturbative aspects of QCD on the CME 
should be quite remarkable. One of the possible non- 
perturbative effects may arise from the instanton-anti- 
instanton (II ) molecules at finite temperature near T c . 

As a major component of the QCD vacuum, the in- 
stantons can produce the quark condensate, the low- 
lying hadron state and other non-perturbative features of 
QCD, according to the interacting instanton liquid model 



(IILM) [11] . At zero and low temperatures, the instanton 
ensemble is in a random liquid state which is responsible 
for the spontaneous chiral symmetry breaking and the 
axial anomaly Naively one would expect that the in- 
stantons might give a little contribution to the CME due 
to their significant suppression at high temperature. 

However, it has been suggested that the hot chiral tran- 
sition is not driven by the vanishing of the instantons 
and anti-instantons but by their rearrangement around 
T c [l2l - [l4j . Namely, when approaching to T c from below, 
the random instantons and anti-instantons are paired up 
into the ordered II molecules. This picture implies that 
the instantons still keeps sizable density near T c with 
the molecular forms. Moreover, a recent investigation 
indicates that the II molecules can even survive up to 

The rearrangement picture indicates that the non- 
perturbative effects induced by the instanton configura- 
tions will persist into the high temperature region. In the 
chiral symmetric phase, the interacting II molecules can 
be regarded as one of the possible mechanisms responsi- 
ble for the sQGP formed in T c < T < 2 - 3T c [lj|. So 
even the chiral restoration is one of the necessary con- 
ditions for the appearance of the electric current j em , 
the instantons may still exert important influence on the 
chiral magnetic effect via the molecular forms. 

It is well known that the random instanton configura- 
tions can lead to the famous 't Hooft interactions among 
quarks with different flavors 

£>tHooft ~ nw°)(^^' w 
/ 

which explicitly breaks the Ua[X) symmetry due to 
the axial anomaly. Similarly, the ordered II molecules 
also result in effective quark interactions, but only 
with the four-fermion coupling forms. For two-flavor 
case, the Fierz-invariant quark interactions induced by 
77 molecules has been derived in [l4[ with the form 

r 2 3 - 

Cmolsym = G[— ^[(V>» 2 - {^^f] 
c a=0 

c c a=0 

+ 0K a 7 „ 75 V0 2 ]} + C 8 , (5) 

where r° and f are unit and Pauli matrices in flavor space 
respectively and G is the coupling constant. The Cg in 
([5]) refers the quark interactions in the color octet chan- 
nel. In the II molecular picture, it is expected that the 
dominant quark interactions for T > T c are the molecular 
ones rather than the 't Hooft interactions [13l - ll6j 2 . Dif- 
ferent from (j4}, the lagrangian ([5]) respects not only the 



2 Note that the 't Hooft type interactions might also survive in 
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SU V (2) <g> SU A (2) ® U B (1) symmetry but also the U A (l) 
symmetry. This is because the II pair gives vanishing 
topological charges [Ill- 
Note that the most general Fierz-invariant form for the 
four-quark interactions under the invariance of SUy (2) C3> 
SUa(2) <8> Ub(1) <& Ua{1) symmetry including both the 
scalar (psudo-scalar) and vector(axial-vector) currents is 

1 3 

- ]r[(^T a 7 ^) 2 + ^T a 7 ^) 2 ] 



a=0 



g^^Y - $rVnsVO a ] + 4 4) , (6) 



where Gi are four independent coupling constants [17| . 
Comparing ([5]) with (|6]), one can find that the interact- 
ing II molecule model (IIMM) gives prediction of these 
coupling constants with 



~ G 3 - J^2 G > ( 7 ) 



which are all dependent on the only parameter G. We 
shall show below that nonzero Gi and G 3 (4) give non- 
trivial contributions to the induced current j em for the 
CME. 

We stress that even we are very interested in the pos- 
sible non-perturbativc influence of the instantons on the 
CME, we do not limit our study within the only for- 
malism ([5]), which is just a special case of ([5]). Without 
loss of generality, we will assume that the Lagrangian 
density with the four-quark interactions (j6|) works for 
T c < T < 2 ~ 3T C . This Lagrangian can be regarded as 
an extension of the traditional NJL model of QCD (In the 
next section, the degree of freedom of the Polyakov-loop 
is also included). 

For the physics related to the CME, we are particu- 
larly interested in the effective four-quark interactions in 
the vector isospin-scalar and axial-vector isospin-scalar 
channels 



Cva = -Gv(ipJ^) 2 - Ga(^7/i75^) S 



(8) 



which appears in both (J5J) and The significance 

of Cva to the CME is attributed to the appearance of 
two special vacuum expectation values (VEVs) in the 
chirally asymmetric system under the influence of an 
external magnetic field B, namely, j z = (ipj 3 ip) and 



the T > T c region according to 1151 . Such interactions should 
have little influence on the CME since the electric current only 
appears in the chirally symmetric phase. 



n,A = (V>7o750}[4]. The former is the VEV of the vec- 
tor current ■07^-0 along the direction of B and the latter 
the chiral charge density conjugate to fj,A- As a conse- 
quence, the chiral magnetic effect might be sensitive to 
the couplings Gy and Ga in (jH]). 

The prediction of the IIMM gives 



G\ 



9i 

2 



G A = 2G 4 + — = -3-^r. (9) 
2N? 2 2N3 



Equation (|9]) indicates two unconventional points in- 
duced by the II molecules: (1) The vector interaction is 
attractive while the negative Ga implies that the axial- 
vector interaction is repulsive; (2) The magnitude of Ga 
is three times that of Gv (For the totally polarized case, 
the coupling Ga in the temporal direction is even 12 
times that ofG v ^). 

This is quite different from the conventional four-quark 
interactions derived from the Fierz transformation of the 
colored quark-antiquark current-current interaction 



Cog 



EA 



(10) 



which arises from the one gluon exchange approximation 
(OGEA) of QCD. In the OGEA, both G v and G A are 
positive with 



r , r G s 
try = <^A = — , 



(11) 



where Gs is the coupling constant in the scalar iso-scalar 
channel. 

Note that it is very likely that the effective quark in- 
teractions stemming from the OGEA can persist into the 
non-perturbative region of QCD with relatively strong 
couplings. A well-known example in the literature is the 
global color model (GCM) of QCD, which is based on 
the nonlocal colored quark-antiquark current-current in- 
teraction 



Cgcm{x,v) = g(x - y)0(x) 7Al A^(a:)0(y)7 M A^(y), 

(12) 

where g(x — y) is a function of the strong coupling [18|. 
This model can successfully describe the low-lying 
hadronspjl and give reasonable QCD condensates in the 
vacuum and at finite temperature and density [20|. 

Since the axial-vector interaction may be repulsive (as 
in the IIMM) or attractive (as in the OGEA or GCM), 
we will treat the coupling Ga in ([8]) as a free parameter. 
The realistic Ga may include contributions from both 
the II molecules and the OGEA, which should depend 
on the temperature in general. On the other hand, since 
the vector interaction is attractive in both the IIMM and 
GCM, we only consider the positive Gv in the following. 
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III. CORRECTION TO THE CME FROM 
AXIAL- VECTOR INTERACTION 



Working at the mean field level, the Lagrangian 
can be rewritten as 



C 



MF 
VA 



G A n A 



2G A n A $ lol5 ib - G v f 2 + 2G v f$j 3 ip. 

We can read two nontrivial points from Eq. (ll3[) . First, 
both squared terms of the condensates j z and n A give 
contributions to the thermodynamical potential. Sec- 
ond, the four-quark interactions in the axial-vector and 
vector channels lead to a dynamical axial chemical po- 
tential \)! A = —2G A n A , and an effective gauge field 
A z = —2Gvj z along the direction of B, respectively. 
Note that both points are not taken into account in [4[ 
to derive the analytical formula for the CME current. 

To investigate the dynamical influence of Cva on the 
CME, we add the kinetic lagrangian for free-interacting 
quarks, namely Ckin — ^{ilfid^ + fi A j°j 5 )tp to the La- 
grangian ^ . We will work in the chiral limit throughout 
this paper. To describe the confinemcnt-dcconfinement 
transition of QCD, we also include the dynamics of the 
Polyakov-Loop, the action of which in the pure gauge 
theory denotes as U. Our model can be looked as the 
Polyakov-loop-enhanced IIMM but with a varying axial- 
vector coupling G a- In form, it is very similar to the so 
called two- flavor PNJL model[2lj|. 

One then gets the mean field thermodynamical poten- 
tial 

ft = U + G s cr 2 + G v f 2 - G A n\ 

AT \if B \ f°° d P* a i \ 

~ n c }^ -k — y, a sk / —f h u s {p) 



f=u 



s,k 



dp z 
2ir 



f=u,d s,k 

x ln(l + 3 $e"^ s + 3 ie^ 2 ^ + e^ 3 ^) , (14) 



where 



S s>+1 for k = 0, eB > , 

<5 s ,'_i for k = 0, eB < , (15) 

l' for k # , 



a s k 



and /a is the UV regulator function. $ and $ in (fT4|) 
correspond to the normalized traced Polyakov loop and 
its Hermitean conjugate respectively. In (fT5"j). s refers to 
the spin and k denotes the Landau level. The dispersion 
relations of the quasi-particles derived from the eigenval- 
ues of the Dirac operator take the form 



where 



uj 2 = M 2 + [|p| + s p. A sgn(p z )]' 



\p z \ 2 + 2\q f B\k 



(16) 



(17) 



and qf stands for the electric charge for /(u or d) quark. 
In the above equations, fi A and p z refer to the effec- 
tive axial chemical potential and the modified momen- 
tum along B respectively, namely 

fj-A = HA + a4 = VA - 2G A n A , (18) 
p z =p z +A Z =p z ~2G v j z . (19) 

The current j z , the chirality density n A , the quark con- 
densate 3 a and the VEV of the Polyakov-loop $ can be 
determined self-consistently by solving the saddle point 
equations 

on an on on 

^r ' ^=°< ^= ' »=°- ^ 

Note that the influence of Gv on the CME has been 
investigated in 22], where only the vector interaction is 
introduced phenomenologically. The main conclusion is 
that Gv may lead to a dielectric correction to the CME 
current density 



1 



3 



2n 2 (l + 2G v C R ) 



N c N f en A B 



(21) 



where Cr — (22[. When including the axial- vector 
interaction, the dispersion relation (fTo| suggests that 
Eq. (f2"Tj) will be further modified as 



1 



J 



2tt 2 (1 + 2G v C r ) 



N c N f eft A B 



(22) 



However, a subtlety on the calculation of the polar- 
ization tensor or susceptibility in [22| has been pointed 
out by Kenji Fukushima (23|: According to Eqs.(8-10) 
in [231 ] . the quantity Cr should be zero since the cor- 
responding static susceptibility II 33 (q° = 0, q — > 0) is 
always vanishing. This implies that the vector interac- 
tion does not give any correction to the CME. Actually, 
this point is also consistent with the Nielsen-Ninomiya's 
argument [24| which is only based on the nonrenormaliza- 
tion of triangle anomalies and the energy conservation. 
As shown in Eq. (l2"Tl) (or Eq.([2"2"|)). the nonzero Cr will 
lead to a modification of the coefficient before \i A B (or 
jl A B in our case). Nevertheless, the argument proposed 
by Nielsen and Ninomiya [24| does not support the co- 
efficient receiving any correction. This is because such a 
coefficient is determined by the the nonrenormalization 
of triangle anomalies, which is known exactly due to the 
topological nature of anomalies (more details of this ar- 
gument see the first part of Sec. Ill in [J] ). 

In light of this, the correct formula for the current den- 
sity should be 



3 = ^N c N f eji A B 



(23) 



3 Strictly speaking, there should be two chiral order parameters in 
the chiral limit since {4>l^Pr) 7^ {^r^l) f° r nonzero fi A . 
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rather than Eg. ([22]) when considering the axial- vector 
interaction. We note that, unlike Eqs. (j2"T1) and (f2"2l . 
Eq. ([23l) does not contradict Nielsen-Ninomiya's argu- 
ment. This is because the energy required to remove 
a particle from the left-handed Fermi-surface and add it 
to the right-handed Fermi-surface is no longer 2/j,a but 
2p>A due to the influence of the axial-vector interaction. 
In fact, Eq. ([2"3")) will be obtained naturally if we take the 
same method as part A of Sec. Ill in [J] but with the 
replacement of Ha with ft a- 

In the following, we will ignore the vector interaction 
and only consider the effect of the axial- vector interaction 
on the CME. Using the method 



Jem 

and the replacement 



an 



A 3 =Q 



d 



dp z 



(24) 



(25) 



one can get the integral expression for the calculation of 
the electric current density in the two-flavor case 



= N C £ 



If 



f=u,d 



kfB \ 

2tt 



s.k 



OCk, 



(26) 

As the case of zero Ga, only the lowest Landau level 
(LLL) gives a non- vanishing contribution in ([2"6)) . Ignor- 
ing the correction from Gy, the induced electric current 
density j em for the two-flavor case under the influence of 
Ga becomes 



N c E 

f=u,d 



2tt 2 



6tt 2 



(27) 



As it should be, Eq. (j2"Tl) indicates that the induced elec- 
tric current density is ultraviolet finite. 

Unlike the analytical formula (|3]), the chirality density 
riA appears directly in Eq. (|2"T1) . This density can be 
obtained self-consistently by 



dn A 



-> riA 



dfJ-A 



djiA 



(28) 



according to (|18l) . This equation is similar to the determi- 
nation of the baryon number density at finite fi when con- 
sidering the vector interaction in the NJL-type model[25[ . 
We note that even the formula (|2"T1) is obtained from the 
LLL integration, it doesn't mean that j em only depends 
on the LLL for nonzero Ga since the high Landau levels 
also give contributions to the chirality density ua- 

The modified expression of j em implies that the in- 
duced electric current density under the influence of the 
axial- vector interaction may deviate from the analytical 
result in the following aspects. First, the current density 
j em is no longer directly proportional to the axial chem- 
ical potential ha and the magnetic field B. The reason 



is that in general, the dynamical chemical potential h'a 
doesn't linearly depend on [iA and B . Second, [i'a is 
T-dependent (Since Ga and ua are both T-dependent), 
so the current density j em also depends on the temper- 
ature. Third, the IIMM predicts that the CME will be 
enhanced significantly: (1) [x' A has the same sign as ha 
because of the negative Ga', (2) h'a may be sizable since 
Ga is strong as Gs- In contrast, the OGEA (or GCM) 
predicts that the CME will be weakened by the attrac- 
tive axial-vector interaction compared to the analytical 
result. 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

In this section, we provide numerical results for the 
influence of the axial-vector interaction on the chirally 
imbalanced matter in a constant magnetic field. The 
coupling Gs, the UV regulator function f\(p) and the 
Polyakov-loop model U in (fl"4f are all adopted from [2(| , 
in which the PNJL model was used to study the CME. 
The model parameters of the NJL part are fitted by re- 
producing the pion decay constant and the quark con- 
densate in the vacuum. The T c for vanishing [iA, B, and 
Ga is ~ 228 MeV with the influence of the Polyakov-loop 
dynamics [26j . 



A. Effect of axial-vector interaction on the chiral 
magnetic effect 

We will focus on the variation of the ratio 



R = jern/jem{GA = 0) = JIa/ Hj 



(29) 



with the impact of the coupling Ga ■ The quantity R = 
R(Ga,T, ha, B) corresponds to the ratio C/2-k 2 , which 
is the focus of the recent lattice QCD calculations [jjllicj. 
The deviation of R from the unity reflects the influence 
of the axial- vector interaction on the CME compared to 
the analytical formula. 

In FigO] we first show the dependence of R on the cou- 
pling Ga for different \iA with fixed eB = 4x 10 18 Gauss 
and T = 230 MeV. The density of II molecules might 
be sizable at this temperature since it is just slightly 
greater than T c (The catalytic effect of the magnetic 
field on the chiral symmetry breaking can be ignored 
for such an eB [26j). We see that R decreases mono- 
tonically with increasing Ga- The electric current den- 
sity is enhanced by Ga < and weakened by Ga > 0. 
Its deviation from the analytical value is quite signifi- 
cant: For the moderate coupling Ga/Gs = —0.5 (0.5), 
the electric current density increases (decreases) by ^62 
(27) percent at fiA = 15 MeV; While for the strong cou- 
pling Ga/Gs — —0.75 (0.75), the electric current density 
increases (decreases) by ^140 (34) percent at the same 
fiA- This is in contrast to the minor modification of j em 
caused by the vector interaction [22j. 
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FIG. 1: The ratio of j em to jem(GA = 0) as a function of Ga 
for different ha- The eB and T are fixed as 4 x 10 18 Gauss 
and 230 MeV, respectively. 



Figure |5] shows the electric current density j em and the 
ratio R as functions of ha for different Ga, where the 
same T and B are used as in FigfT] For Ga = —0.5, 
the upper panel indicates that the electric current den- 
sity is nearly proportional to ha in the small ha region, 
but its deviation from the linear behavior becomes more 
and more significantly with increasing ha ■ Different from 
this, the electric current density almost linearly depends 
on ha up to the moderate ha region for Ga = 0.5. These 
features are demonstrated more clearly in the lower panel 
of Figf2] The ratio R monotonically increases with ha 
for Ga = —0.5 while almost remains as a constant for 
Ga = 0.5. The same information can be read from Fig|T] 
This is because the chirality density rises more rapidly 
with growing ha for the repulsive axial-vector interac- 
tion. 

The dependence of R on the temperature with ha = 
15 MeV and eB = 2x 10 18 Gauss is displayed in FigH 
For Ga < 0, the ratio R monotonically increases with 
increasing T and its deviation from one becomes huge for 
the strong coupling. For Ga > 0, the ratio R decreases 
relatively slowly with growing T for both the moderate 
and strong couplings. Though the change of R by T is 
not so significant as Ga < 0, the ratio R still decreases by 
-10 (15) percent from ~ T c to ~ 1.1T C for G A /G S = 0.5 
(1). Hence, the temperature has a noticeable effect on 
the induced current. The reason is that the density ha 
is a monotonically increasing function of T, as reported 
in HI. 

Figure|4] indicates that the ratio R is not so sen- 
sitive to the external magnetic field but still slowly 
rises(drops) with eB for Ga < 0(Ga > 0). This is be- 
cause that the chirality density is enhanced by the mag- 
netic field through the increase of the thermodynamical 
potential^- Note that such a tendency will slow down 
or cease for large enough B since its catalytic effect on 
the chiral symmetry breaking: the enhancement of the 
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FIG. 2: Upper panel : The current j em as a function of ha for 
different Ga- Lower panel: The ratio of jem to j e m{GA = 0) 
as a function of ha for different Ga- The eB and T are fixed 
as 4 x 10 18 Gauss and 230 MeV, respectively. 



chiral condensate by the magnetic field will suppress the 
chirality density. 

Figures. 1-4 imply that the chirality density-density 
correlation has an important impact on the chiral mag- 
netic effect, no matter the axial- vector interaction arises 
from the 77 pairings, the OGEA, or the GCM. The ra- 
tio R can deviate significantly from the unity under the 
combined effects of Ga, Ha, T and B. The deviation may 
be positive or negative which depends on the sign of Ga- 
This contrasts with the role of Gy on the CME: It only 
affects the CME obviously when the magnetic field B is 
very large [12 i Its correction to the CME has nothing to 
do with ha and T. 

Let us compare our results with the recent lattice QCD 
data. In @, the CME was qualitatively confirmed by the 
lattice simulation via introducing a finite ha- However, 
the obtained coefficient C is around 0.013, which is far 
less than the analytical value l/27r 2 ~0.05. The further 
study using the quenched simulation suggests that the 
suppression is sensitive to the lattice spacing[10]. The 
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for different Ga- The ha and T are fixed as 100 MeV and 
240 MeV, respectively. 



new lattice coefficient C is roughly in the range 0.02-0.03 
10], which increases greatly compared to the former re- 
sults. However, it still much less than the standard ana- 
lytical value. Moreover, another new point in [l(| is that 
the coefficient C might be dependent on the temperature. 

Regardless the possible errors in .10] , we stress that the 
current lattice data are quite natural when considering 
the effect of the attractive axial- vector interaction: (1) 
According to (|27|l , the induced electric current is propor- 
tional to the effective axial chemical potential p,Ai which 
may be suppressed significantly compared to fiA by the 
positive Ga- As demonstrated in Figs. 1-4, the suppres- 
sion is really remarkable for the moderate and strong 
couplings. So the reduction of \1a by the attractive axial- 
vector interaction can explain why the lattice data of the 



induced electric current are much less than the analytical 
value. Let us make a rough comparison. According to 
Fig.7 in pi], the lattice ratio C/2tt 2 is located in 0.4-0.6 
for T ~ 1.1T C (T C ~ 270 MeV). Our numerical results 
shown in Fig J3] suggest that the ratio R is in the range of 
0.5 - 0.65 for T ~ 1.1T C (T C - 230 MeV) if G A /G S varies 
from 1 to 0.5. We see that the lattice data are quite 
consistent with our results if the axial-vector interaction 
is attractive and strong enough. (2) The lattice data 
suggest that the coefficient C or the ratio C/2tt 2 is in- 
sensitive to fj,A &ud eB. This feature is also in agreement 
with our results for the attractive axial- vector interac- 
tion, as shown in Figj2] and Figj4] (3) One set of data 
{/3 = 5.90 with N t = 4, 6, 12) in Fig.7 of [13] indicates 
that the coefficient C decreases with growing T. The au- 
thor has argued that such a behavior might be a lattice 
artifact. However, that the induced current density j em 
may be suppressed by T can be well interpreted by the 
attractive coupling Ga as demonstrated in Fig(3J Ac- 
cording to Fig.7 of . the ratio C/2ir 2 decreases from 
~ 0.53 to ~ 0.32 when increasing T from 317 MeV to 
475 MeV; Or in other words, the ratio C/2n 2 decreases 
roughly by ~ 0.04 per 0.1T C . This is not so far away from 
the decline rate ~ 0.07 per 0.1T C for Ga/Gs — 0.5, as 
indicated in Fig(3] 4 So to judge whether the coefficient C 
is sensitive to the temperature or not needs more detailed 
lattice investigation. 



On the contrary, the IIMM predicts that the CME 
current or the coefficient C will be enhanced by the 
II molecules. Clearly, this prediction deviates distinctly 
from the present lattice data. This might suggest that the 
dominant axial-vector interaction for T > T c is coming 
from the OGEA. Nevertheless, the current lattice simu- 
lation is still rough and far from conclusive(at leat at the 
quantitative level). First, that the current lattice data 
are sensitive to the lattice spacing [l(| suggests that the 
continuum limit is very important for quantitatively un- 
derstanding the chiral magnetic effect. Second, the main 
simulations in jlQ] are still at the quenched level while 
the discretization error of the fermion action is larger 
than that of the gauge action. Therefore, more improved 
quantitative lattice investigation will shed light on the 
role of the II molecules on the chiral magnetic effect. 



In any case, we stress that, due to the contribution of 
the dynamical axial chemical potential, the lattice result 
of the induced current for a finite \xa doesn't need to be 
quantitatively consistent with the analytical formula. If 
no other QCD modifications, the coefficient C or the ratio 
R should be greater (less) than the standard analytical 
value for G A <0 (G A > 0). 



4 In the realistic case, the coefficient C should change more mildly 
since the coupling Ga becomes weaker with increasing T. 
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FIG. 5: The locations of the chiral TCP in the T — /ia plane 
for different Ga- 



B. Effect of axial- vector interaction on the T — ha 
phase diagram 

The chiral tricritical point (TCP) on the T — /j,a phase 
diagram has been found in [26| , which locates at the rel- 
atively larger [ia region (in the chiral limit for vanishing 
quark number density). The TCP is also confirmed in 
the linear sigma model, whose location is very close to 
the T-axis|27[. These studies show that the T — fi A phase 
diagram is somehow similar to the T — fi phase diagram 
of QCD. 

Here we will demonstrate the role of the axial- vector 
interaction on the T — fiA phase diagram. We are partic- 
ularly interested in its influence on the TCP. In Fig|5j we 
show the locations of the TCP for different axial- vector 
interactions. We see that the TCP is quite sensitive to 
the axial- vector coupling. The negative Ga shifts the 
TCP towards the smaller [ia region while the positive Ga 
moves it towards the larger fiA area. This implies that 
the first-order chiral transition is strengthened by the re- 
pulsive axial- vector interaction induced by // molecules 
while weakened by the attractive axial- vector interaction 
stemming from the OGEA. 

It is well-known that the attractive vector interaction 
can soften the chiral transition and move the critical 
point(CP) towards the larger [i area via generating a dy- 
namical quark chemical potential [111 [28| . Especially, 
when taking into account the two-flavor color super- 
conductivity, the vector interaction combined with the 
charge- neutrality constraint and/or the axial anomaly 
can even lead to multiple chiral critical points 29] 5 . We 



5 The new critical point might also appear through the inter- 
play between the chiral condensates and the diquark conden- 
sates in the color-flavor-locking phase in the context of the axial 
anomaly ;30|. 



observe from Fig JS] that the attractive axial- vector inter- 
action at finite fiA has a similar influence on the chiral 
phase transition as the vector interaction does at finite 
fi. 

However, in the instanton molecule model, the axial- 
vector interaction at finite fiA and the vector interaction 
at finite fi have opposite effects on the chiral transition. 
Since the induced Ga and Gv have different signs, the 
chiral transition will be strengthened (softened) at finite 
fj,A (p). In addition, the influence of the axial- vector 
interaction on the chiral transition at finite \ia could be 
more remarkable since the coupling G a is much stronger 
than Gy in the instanton molecule picture. 

Locating the CP on the T — fi phase diagram is a fo- 
cus of the current studies on the QCD phase transition 
and heavy-ion collisions. So far, the location and even 
the existence of the CP is still under debate. Since there 
is no sign problem at finite [ia, it is very interesting to 
explore the T — \ia phase diagram of the chirally imbal- 
anced matter and locate the TCP by means of the lattice 
simulation. Such an investigation may provide useful in- 
formation on the true QCD phase transition. Moreover, 
studying the T — /ia phase diagram through lattice QCD 
may shed light on the non-perturbative features of sQGP. 
For example, we can test the instanton molecule picture 
of QCD by investigating the T — /j,a phase diagram. 



V. SUMMARY AND CONCLUSION 

In this article, we have first investigated the influence 
of the axial- vector interaction on the chiral magnetic ef- 
fect by introducing an axial chemical potential ha- Such 
an interaction can be induced by the instanton molecules 
or derived from the one gluon exchange approximation 
of QCD. In the presence of a finite chirality density ha, 
a dynamical axial chemical potential fj,' A is generated 
through the density-density correlation. We derived that 
the induced electric current density j em in an external 
magnetic field B linearly depends on the effective axial 
chemical potential \1a which including the contribution 
of \J A . Accordingly, the obtained j em deviates from the 
standard analytical formula and explicitly depends on 
the chirality density ha- In general, the deviation relies 
on the axial- vector coupling Ga, the axial chemical po- 
tential (ma, the temperature T, and the magnetic field 
B, 

For Ga > from the one gluon exchange approxi- 
mation, the induced current density can be significantly 
suppressed comparing with the analytical value. We also 
find that the ratio R — jem/ jem(GA — 0) is insensitive to 
the axial chemical potential and the magnetic field but 
decreases with increasing temperature if Ga > 0. These 
features are quite consistent with the recent lattice results 
at finite [ia- Hence, the suppression of the lattice data 
comparing with the analytical formula can be attributed 
to the influence of the attractive axial-vector interaction 
of QCD. Actually, it is very likely that the deviation of 



9 



the lattice result originates from the difference between 
Ha and ha- The reason is that Nielsen- Ninomiya's ar- 
gument does not support the renormalization of CME if 
HA is not shifted dynamically. 

On the contrary, an unconventional prediction from 
the instanton molecule model is that the axial-vector in- 
teraction is repulsive with G a < and much stronger 
than the vector interaction. As a consequence, the chi- 
ral magnetic effect will be enhanced significantly by the 
instanton molecules near T c . In such a picture, the ratio 
Jem 1 1 jem{G ' a = 0) increases with both \xa and T. These 
features deviate obviously from the current lattice re- 
sults. Since the present lattice simulation at finite ha is 
still rough and sensitive to the lattice spacing, we antic- 
ipate the more improved lattice calculation to test these 
predictions. 

In addition, we also demonstrated that the axial- vector 
interaction plays an important role on the T — ha phase 
diagram of the chirally imbalanced matter. The repul- 
sive (attractive) axial-vector interaction can effectively 
strengthen (weaken) the chiral phase transition if the 
coupling is relatively strong. The chiral TCP on the 
T — [i A plane is shifted towards the lower (higher) ha re- 
gion by the repulsive (attractive) axial-vector coupling. 
This is quite different from (similar to) the role of the 
vector interaction on the chiral phase transition at finite 
baryon chemical potential /i. 

Since there is no sign problem at finite ha, the im- 
proved lattice simulation is still very welcome for dis- 
closing the non-perturbative QCD effect on the chiral 
magnetic effect. Moreover, studying the deviation of the 
induced current density from the analytic result and lo- 
cating the TCP in T — ha plane by the lattice simulation 
can be used to test the instanton molecule picture for 
T>T C . 

We note that if ha is introduced by coupling it only to 
the chiral density, it can't be considered as a true chem- 
ical potential. The reason is that the axial charge is not 
a conserved quantity due to the axial anomaly. Recently, 
it is proposed by Robakov [3l| that ha should be conju- 



gated to a proper combination of the chiral density and 
a Chcrn-Simon term, which is a conserved quantity. It is 
then claimed in [3l[ that the fact ha must be associated 
with a conserved charge is essential in the discussion for 
the nonrenormalization of the CME current. 

Since there exists a subtlety on the definition of ha, 
a question naturally arises that whether the dynamical 
axial chemical potential h'a defined in this paper can be 
associated with a conserved charge. We point out that 
h'a doesn't need to be conjugated to a conserved charge if 
Ha could be really well-defined. The reason is that h'a ^ s 
an induced quantity generated dynamically by the axial- 
vector interaction for nonzero ua- This is different from 
the parameter ha (Assuming it has been well-defined as a 
true chemical potential) , which should be associated with 
a conserved charge in principle. We note that no matter 
how ha is introduced, it must give rise to nonzero ha- 
As demonstrated in this paper, the h'a induced by ua 
will modify the chiral imbalance described originally by 
HA- 

A natural extension of this work is to investigate the 
role of the axial- vector interaction at finite h i n the pres- 
ence of an external magnetic field. Conjugate to the 
induced current for the chiral magnetic effect at finite 
Ha i an axial- vector current along the direction of the ex- 
ternal magnetic field is generated at finite ^[32||. An- 
other related topic at finite /i is the so called chiral 
shift parameter 33] in a magnetic field. The roles of the 
axial-vector interaction in these cases will be reported 
elsewhere [13] ■ 
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